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ABSTRACT 

We report observations of the dust cloud L1204 with the Onsala 25-m tele- 
scope in the 6 cm (I11-I10) transition of H 2 CO. The observed region includes the 
S140 Ha arc. This spectral line is seen here in absorption against the cosmic 
microwave background, indicating the presence of widespread warm molecular 
gas at intermediate densities. Overall, the distributions of H 2 CO and CO (taken 
from the literature) are fairly similar, though significant differences exist at small 
scales. Most notably, while the CO peak is nearly coincident with the S140 Ha 
arc, the maximum H2CO absorption is clearly separated from it by a full 10' beam 
(~ 3 pc). We argue that these differences result from differing abundances and ex- 
citation requirements. The CO(l-O) line is more optically thick and more biased 
towards warm gas than the H2CO 6 cm line. On the other hand, formaldehyde is 
more easily photodissociated and is, therefore, a poorer tracer of the molecular 
gas located immediately behind Photon Dominated Regions. 



Subject headings: ISM: clouds — ISM: molecules - radio lines: ISM - stars: 
formation - galaxies: ISM 
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1. Introduction 

Since H2 -the most abundant molecule in space- lacks a permanent dipole moment, its 
rotational transitions are prohibited. Although the quadrupolar transitions exist, they are 
of little use for the syudy of the bulk of molecular gas in the ISM because they require high 
temperature to be excited. Instead, the structure and properties of cold molecular clouds 
in the interstellar medium are usually studied using low-energy rotational transitions of 
simple non-symmetric polar molecules. For practical reasons, the first rotational transition 
(J = 1— »0) of carbon monoxide (CO), at 115.27 GHz has been the most popular choice. 
This transition, however, has long been known to be nearly always optically thick, so - 
for a given filling factor- its intensity is expected to increase monotonically with the kinetic 
temperature of the emitting gas. Clearly, this could have adverse effects on efforts to establish 
the distribution of molecular gas in the ISM from CO observations alone, because very 
low temperature gas might go unnoticed in sensitivity-limited CO observations while warm 
regions (> 20 K) will stand out even if they are not those with the highest molecular content. 
The sources where these effects might be most noticeable are those with large temperature 
gradients; for instance in molecular clouds located in the immediate vicinity of hot stars. 

While all molecular emission tracers share this temperature dependance to some degree, 
absorption lines can be detected even in very cold gas, provided sufficiently bright background 
continuum sources are available. The scarcity of such sources at the wavelength of the 
common molecular tracers, however, has limited the usefulness of absorption measurements 
in the study of specific Galactic molecular clouds (e.g Evans et al. 1980). The 6-cm (I10-I11) 
transition of ortho formaldehyde (H 2 CO) offers an interesting alternative. Owing to collisions 
with neutral particles that selectively overpopulate the lo wer energy level, the exc itation 



temperature of the li — > In transition lies below 2.7 K (ITownes fc Cheung! 119691 ). This 



allows t he transition to be o bserved in absorption against the cosmic microwave background 



(CMB) (jSnyder et al.lll969l ). and makes it a potentially powerful tracer of molecular gas in 
any direction of the sky. The excitation requirements are such the the 6 cm H 2 CO line is a 
good indicator of the presence of cool to warm molecular gas (T > 10 K) at intermediate 
densities (10 3 cm~ 3 < n < 10 5 cm -3 ). Unfortunately, the absorption line is weak, so very 
large amounts of telescope time are required to map large areas of the sky. 



Recently, iRodrfguez et al.l (120061 ) conducted a blind search for H 2 CO absorption and 
compared CO emission and H 2 CO absorption profiles towards the Galactic anticenter. They 
found a rough, large-scale correlation between these two tracers, and concluded that both 
lines preferentially trace warm and dense molecular gas. Here, we will examine this relation 
between CO and H 2 CO at a somewhat smaller scale using observations of the well-known, 
nearby star-forming region Sharpless 140 (S140 -Sharpless 1959) associated with the dark 
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dust cloud Lynds 1204 (L1204 -Lynds 1962). L1204 is centered at I = 107?47, b = +4?82 
and covers an area of 2.5 square degrees (Lynds 1962). At its southwest edge lies S140, a 
prominent compact arc-shaped Hll region with an angul ar size of ~ 2' x 6 '. The ionization 
of S140 is maintained by the nearby BOV star HD211880 fakir et ali |l9 7_8]). The distance of 
S140 /L1204 deduced from the brightness of the exciting star is 910 pc (ICrampton fc Fisher 



1974 ). S140 has been the subject of many observational studies, that have usually focused on 



the Photon Dominated Region (PDR) on the edge of L1204, and on the embedded infrared 
sources located right behind it (e.g., Preibisch et al. 2001; Hayashi & Murata 1992; Preibisch 
& Smith 2002; Bally et al. 2002). Remarkably, while the dust cloud is seen as an extended 
dark feature covering more than two square degrees, the CO emission peaks immediately 
behind the Ha arc (Heyer et al. 1996; Evans et al. 1987; Blair et al. 1978), while only 
relatively faint emission extends deep within the dust cloud (jHelfer fc Blitzlll997l ) . The 6-cm 



l ine of H 2 CO was detected in absorption against the CMB i n L1204 near S140 b y lBlair et al. 



(119781 ) with the NRAO 43 m telescope, and unexpectedly by Evans et al.l (119871 ) during VLA 
observations of the bright condensation just northwest of S140. However neither of those 
studies provided a extensive mapping of the H 2 CO CMB absorption in L1204, and the exact 
extension of the gas traced by H 2 CO remains unclear. In this article, we will present such a 
extensive mapping of the 6 cm CMB absorption of H 2 CO over most of the large dust complex 
L1204, and compare our results with existing CO observations taken from the literature. 



2. Data 

The H 2 CO observations were obtained during two sessions (January and September- 
October 2004, respectively) with the 25.6-m telescope of the Onsala Space Observatory 
(OSO) in Sweden. At 6 cm, the angular resolution is 10', and our pointing precision was 
always better than 20". Frequency-switching, with a frequency throw of 0.4 MHz was used, 
and both polarizations of the incoming signal were recorded simultaneously in two indepen- 
dent units of the autocorrelation spectrometer. Each of these units provided 800 2 kHz-wide 
channels. At the observed frequency of 4829.660 MHz, this setup provided a total bandwidth 
of 99 km s _1 and a (Hanning-smoothed) velocity resolution of 8 kHz = 0.49 km s . The 
spectrometer was centered at the systemic velocity of S140, Vlsr = —8.0 km s _1 . Daily 
observations of the supernova remnant Cas A were used to check the overall performance of 
the system. The system temperature during our observations varied from 33 to 36 K. 

In order to map the entire region behind S140, we observed 72 positions on a regular 
square grid with a 10' spacing, centered at I = 107P0, b = +5?3; the resulting map uniformly 
covers a 1?0 x 1?8 rectangular region (Fig. 1) . The off-line data reduction was done 
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with the CLASS program of the GILDAS software package (iGuilloteau fc Forveilld Il989l ). 
and involved only the subtraction of (flat) baselines from individual integrations and the 
averaging of all spectra taken at the same pointing position. The total integration time for 
each of these positions was about 10 hours, yielding a typical final noise level of 3 mK (T^). 

The distribution of radio continuum sources in the region of L1204 has been studied in 
detail by Allen Machalek & Jia (in preparation), using data from the Canadian Galactic Plane 
Survey. Fairly bright continuum emission is associated with the Ha arc and the embedded 
massive protostars located behind it -but, as we will see momentarily no formaldehyde 
was detected from either of these regions. In addition, a number of extragalactic background 
sources as well as diffuse emission associated with the dust cloud L1204 itself contribute to the 
overall radio continuum. The typical brightness temperature average over the Onsala beam 
at 6 cm, however, is only about 0.2 K, except towards the Ha arc and the embedded massive 
protostars (where again, no absorption was detected). Since the brightness temperature is 
so small, any H 2 CO absorption profiles features detected must be absorption of the cosmic 
microwave background radiation at 2.7 K. 

In the analysis of our new observations, we will also make use of 12 CO(1-0) observations 
of L1204/S140 kindly provided by Dr. Tamara Heifer, and published in Heyer et al. (1996), 
and Heifer & Blitz (1997). These data were obtained with the 14-m telescope of the Five 
College Radio Astronomy Observatory (FCRAO) in Amherst (MA), and have an intrinsic 
angular resolution of 45". For comparison with our formaldehyde data, we have smoothed 
the CO (1-0) observations to 10', and resampled them on our observing grid. 



3. Results 

Formaldehyde absorption was detected in at least 16 of our 72 observed positions (see 
Fig. [2]). The maximum absorption is located 10' arcmin behind the S140 Hll region at a 
LSR velocity of —8.0 km s -1 , similar to that of the CO emission detected in that area. 



Table 1. Source positions. 



Source 


Position (/, b) 


Size (Aa, AS) 


a(J2000.0), 5(J2000.0) 


Reference 


L1204 


107?37, +4?87 


L°0 x 


2?5 


22 h 26 m 4i. s +63 15'36" 


Lynds (1962) 


S140(H a ) 


106°8, +5?3 


2' x 


6' 


22 h 19 m 23. s +63°18'16" 


Sharpless (1959) 


Our survey 


107P0, +5?3 


lfO x 


1?8 


22 h 20 m 52. s +63°24'49" 


This paper 
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Right Ascension 
22 b :24 m :0Q 5 22:22 22:20 22:18 22:18 22:16 22:14 




Galactic Longitude 



Fig. 1. — This figure shows a sketch of the L1204/S140 region, with our observed positions 
shown as " +" symbols. The correspondence between Galactic coordinates (used throughout 
the paper) and equatorial coordinates (that have usually been preferred for observations of 
S140) is indicated. The circle represents the 2?5 size of the dust cloud L1204 as reported 
by Lynds (1962). The central positions of L1204 and S140 are shown as "" symbols. The 
contour represents the lowest value of the H 2 CO absorption. 



A second spatio-kinematical structure is detected towards the north-east (here, and in the 
rest of the paper, north and all other directions refer to Galactic coordinates), at Vlsr ~ 
— 11 km s -1 . Both components are presumably associated with L1204, and have clear CO 
counterparts (Fig. [2]- Blair et al. 1978; Evans et al. 1987; Sugitani & Fukui 1987; Park & 
Minh 1995). There is also an isolated absorption feature towards the southeast, at Vlsr ~ 
—2.5 km s -1 . Given its low LSR velocity, this feature i s likely unrelated t o L12 04, and is 
probably a local cloud along the line of sight. Thus, while ISugitani fc Fukuil (119871 ) identified 
three molecular components associated with L1204 in their 13 CO observations, we only find 
two in our formaldehyde d ata. We do find ev idence, however, for a systematic velocity 
gradient across the cloud. iPark fc Minhl (119951 ) argued that this complex overall spatio- 



kinematical morphology was created when S140 and L1204 were swept up by an expanding 
shell associated with the Cepheus bubble. Our data do not illuminate this assertion any 
further, and a more thorough study is necessary to understand the detailed structure of this 
region. 
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Galatic Longitude Offset (arcir.in) 



Galactic Longitude Offset (aremm) 



Fig. 2. — (a) Mosaic of H2CO CMB absorption spectra observed in the L1204/S140 region. 
The (0,0) position corresponds to I = 107?0, b = 5?3, about 12' east of S140, and about 
40' north-west of the nominal center of L1204 (Lynds 1962). (b) Corresponding CO(l-O) 
observations smoothed to 10' (see text). Note that the CO emission is located a full 10' west 
of the maximum H 2 CO absorption. 



4. Comparison with other molecular tracers 



The S140/L1204 region has been observed in many different molecular tracers (e.g. 
Tafalla et al. 1993, Zhou et al. 1993, Park & Minh 1995), but most of these observations 
have focused either on the S140 PDR or on the embedded infrared sources located just 
behind S140, while only a few observations covered the entire dust cloud. Indeed, the first 
CO observations of S140 (iBlair et aJj Il978l ) only covered a limited part of the region. To 
our knowledge, the only existing large-scale CO map of L1204 is that obtained in the 90s 
with the FCRAO telescope (see §2) and published by Heyer et al. (1996) and Heifer & Blitz 
(1997)0. As mentioned earlier, we will use a smoothed version of that dataset here in order 
to compare with our formaldehyde observations. 

In general, the CO emission and H 2 CO abs orption mor p holog ies in this region are quite 
similar (Fig. [3]). This was already noticed by iBlair et al.l (119781 ) in their 6' observations. 
It is also in good agreement with the results obtained towards the Galactic anticenter by 



1 The region lies on the edge of, and is only partly covered by, the CfA Galactic plane survey of Dame et 
al. (2001). 



- 7- 



H,CO 
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107.6 107.4 107.2 107 106.8 106.6 107.6 107.4 107.2 107 106.8 106.6 
Galactic Longitude Galactic Longitude 



Fig. 3. — All panels show a grey-scale version of the DSS-red image of the region around 
/ = 107°0, b = 5f4. The H a arc of S140 is clearly visible on these (red) images near I = 
106P8, b = 5? 3. In the thr ee left panels (al-a3), CO contours taken from the smoothed CO 
data of iHeyer et al.l (119961 ) are overlaid on top of the DSS image, whereas in the three right 
panels (bl-b3), our H 2 CO contours are overlaid. The contours in the top two panels (al-bl) 
include the entire velocity range associated with L1204 (from -12 to -5 K km s _1 ), whereas 
in the middle two (a2-b2) and bottom two (a3-b3) panels, the contours correspond only to 
the -11 K km s^ 1 and the -8 K km s~ x components, respectively. The asterisks correspond 
to the position of the peak H 2 CO absorption for each component. 
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Rodrfguez et al. (2006), and towards the Orion molecular complex by Cohen et al. (1983). 
There are, however, several noteworthy differences between the CO emission and H 2 CO 
absorption in S140. The first difference is the fact that the CO peak and the H 2 CO absorption 
maximum are not located at the same position. The CO integrated intensity map (Fig. [31 
see also Fig. [2]) shows that the maximum CO emission occurs just behind the S140 Ha arc 
at the western edge of L1204, while only comparatively fainter emission extends to greater 
longitudes. The maximum H 2 CO absorption, however, is located about 10' eastward of the 
CO peak. A second notable difference between the CO and H 2 CO profile is the existence of 
a CO emission "tail" in the south/southeast part of the main cloud with little or no H 2 CO 
counterpart (Figs. 2 and 3). 



Intensity ratio 



O 
o 



S140 + L1204 
L1204 

Galactic Anticenter 




I ( CO) K km/sec 



Fig. 4. — Correlation between the H 2 CO absorption line intensity and the 12 CO(1-0) emission 
line intensity at corresponding points in L1204. The squares correspond to data from Table 
El the open circles correspond to H 2 CO upper li mits and the triangles correspond to the 
Galactic anticenter data previously published in iRodrfguez et al.l (120061 ) . The horizontal 
and vertical lines show the "best case" detection limit. The dashed line is the least-squares 
fit for the L1204/S140 region data, and the dash-dotted line is the least-squares fit for the 
Galactic anticenter data from Rodriguez et al. (2006). Note that the fits do not pass through 
the (0,0) point, suggesting that the relation is not linear at low intensity values. 



In order to study the relation between H 2 CO absorption and CO (1-0) emission in a 
more quantitative way, we have computed the moments of the profiles shown in Fig. [2J The 
results are listed in Table [2] (Appendix B). When two velocity components are visible at 
a given pointing, the moments for each were computed separately. Intensities above 3a are 
shown as squares in Fig. [U and were used to make least-square fits (see below). The two 
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spatio-kinematical components that we identified in our formaldehyde dataset behave quite 
similarly with respect to the CO-H 2 CO relation, and are plotted together in Fig. HI The 
best least-squares fit to a straight line for the entire L1204 dataset yields: 

I(H 2 CO) = (4.4 ± 1.1) x 1(T 3 I (CO) + (34 ± 16) x KT 3 K km s" 1 . (1) 

We shall see momentarily that the CO emission near S140 may be particularly bright because 
of local heating. Ignoring the pointings very near S140, however, yields a fairly similar 
relation between CO and H2CO: 

I(H 2 CO) = (3.8 ± 0.9) x 1(T 3 I(CO) + (31 ± 9) x 1(T 3 K km s" 1 . (2) 

The small difference between these relations presumably reflects the differing excitation 
requirements for the two lines. The relation between CO and H2CO given by Eqs. [1] and 
[2] for the L1204 region is almost identical to that found towards the Galactic anticenter by 
Rodriguez et al. (2006): 

I(H 2 CO) = (4.1 ± 0.5) x 10~ 3 I(CO) + (34 ± 13) x 10" 3 K km s" 1 . (3) 

It is important to note, however, that, in spite of the agreement between the fits to 
the Galactic anticenter and S140 data, there is very significant scatter in the CO-H 2 CO 
relation, some poi nts lying nearly lOcr aw ay from the linear relation. This situation was 



already noticed by IRodriguez et all (120061 ) in their study of the Galactic anticenter. This 
lack of a detailed correspondence between CO emission and H2CO absorption presumably 
reflects differences in the excitation conditions of the two tracers, as we will now discuss in 
the next section. 



5. Discussion 

The comparison between H2CO absorption and CO(1-0) emission profiles in the Galac- 
tic Anticenter and in the L1204/S140 region has led us to three important observational 
conclusions: 



1. Qualitatively, the morphology of CO and H 2 CO are quite similar, and quantitatively, 
the line integrated intensities correlate quite well with one another. 
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2. The scatter in the CO-H2CO relation is, however, significantly larger than the obser- 
vational errors. 

3. In the specific case of S140, the CO emission peak is offset by about 3 pc from the 
locus of the deepest formaldehyde absorption, and there is a region south of the main 
cloud where significant CO emission is detected with little or no H 2 CO counterpart. 



From the general large-scale correspondance between the CO (1-0) and H2CO 6 cm 
integrated maps (Fig. [3]), and from the fair correlation between their line intensities (Fig. 
HJ), we conclude that the physical conditions needed for the excitation of both lines are 
quite similar. The calculations presented in the Appendix A, indeed show that both lines 
preferentially trace warm gas at intermediate densities (10 3 6 < n < 10 5 for H 2 CO; n > 10 2 ' 8 
for CO). In this scheme, the offset between the CO(l-O) peak and the H 2 CO maximum 
absorption may seem puzzling. Note that a similar trend is seen at higher resolution: while 
the CO peak in the full-resolution CO map published by Heyer et al. (1996) and Heifer & 
Blitz (1997) is at I = 106°. 8, b = +5°. 3, the H 2 CO absorption feature seen in the high- 
resolution VLA images published by Evans et al. (1987) is centered around I = 106°. 9, b = 
+5°. 3, again a few arcminutes to the east. 

We suggest that a combination of two effects may explain this puzzling result. First, it 
can be seen from the excitation analysis presented in Appendix A that the H 2 CO absorption 
strength "saturates" for T > 30 K, whereas the temperature of the CO emission continues 
to rise at higher kinetic temperature (Fig. 5). For example, while the H 2 CO line strength 
increases by only about 30% when the kinetic temperature goes fro m 20 to 40 K, the CO (l-O) 



line intensity increases by more than a factor of two. According to lPark fc Minhl (119951 ). the 
CO brightness temperature is about 40 K at the peak and 20 K for the rest of the cloud. Thus, 
the strong CO peak behind S140 may well be largely due to enhanced kinetic temperatures 
related to local heating (by the external star providing the ionization of S140, and/or by 
the infrared sources embedded in the cloud). As one progresses into the cloud, the local 
heating diminishes, and the CO line intensity fades. Also, it should be pointed out that the 
formaldehyde calculations presented in the Appendix show that the 6-cm line should be seen 
in emission rather than absorption when the density exceeds 10 5 cm -3 . The fact that this is 
not the case near the CO peak (neither in our low-resolution data, nor in the high-resolution 
VLA data presented by Evans 1978) suggests that the gas density there is lower than 10 5 
cm -3 . Other effects that could explain the offset between the CO and the H 2 CO peaks are 
the lower dissociation energy, and the lower abundance (and, therefore, lower self-shielding) 
of formaldehyde compared to CO (see Appendix A. 3). In a photo-dissociated region, these 
effects should combine to create a stratified distribution where CO survives nearer the source 
of the UV photons than H 2 CO. This stratification, combined with the heating of the CO, 
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would naturally lead to the offset between H 2 CO and CO seen in the present data. 

Finally, the origin of the other main difference between CO and H 2 CO in S140, namely 
the existence of CO emission at the south of L1204 with no or little formaldehyde counterpart, 
is likely related to another aspect of the excitation differences between the 6 cm line of 
formaldehyde and the 1-0 transition of carbon monoxide. Fig. 6 of the Appendix A. 2 shows 
that the density detection limit for H 2 CO line is ~10 times larger than the density limit for 
the CO (1-0) line. We therefore suggest that the gas traced by the CO emission to the south 
of L1204 is of relatively very low density. It is interesting to note, indeed, that classical 
high-density molecular tracers (e.g. CS or NH3) have only been detected around the CO 
peak behind S140, and not in the southern region of the cloud. 

Thus, we conclude that the CO(l-O) and H 2 CO 6 cm lines both tend to preferentially 
trace warm gas at intermediate densities. There are, however, significant differences related 
either to differing excitation requirements or to differing abundances. These differences can 
easily explain the large scatter in the CO-H 2 CO relation. 

6. Conclusions 

The main conclusions of this work are the following: 

1. We have mapped a large region (70' x 110') around L1204/S140 in the 6 cm line of 
formaldehyde, observing a total of 72 regularly-spaced positions every 10' on a regular 
grid. The center of our map was at / = 107?0, b = +5?3, and formaldehyde was 
detected against the cosmic microwave background in at least 16 of our 72 positions 
(Fig. 2). 

2. The formaldehyde emission can be separated in three spatio-kinematical components 
(Fig. 3): two (at Vlsr ~ -11 km s _1 in the northeast part of the cloud, and at Vlsr 
~ -8 km s~ x just behind S140) are clearly associated with L1204, whereas the other 
(an isolated component at Vlsr ~ ~2.5 km s" 1 towards the southeast) is most likely a 
local foreground cloud unrelated to S140/L1204. 

3. Both qualitatively and quantitatively, the CO(1-0) emission and the formaldehyde 
6 cm absorption lines correlate fairly well. An excitation analysis shows that both 
preferentially trace warm gas at intermediate densities. 

4. There are, however, notable differences between the CO and H 2 CO lines, that can be 
traced to differing excitation requirements and abundances. Those differences are most 
likely the origin of the large scatter in the CO-H 2 CO intensity correlation. 
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Model calculations 



A.l. Collisional pumping 



Observations of H2CO in dark clouds show that the anomalous absorption of the 6 and 
2 cm lines is due to col lisions with H2 that selectively overpopulate the lower levels of the 
lines (jEvans et al.lll975l ). We calculated the non-LTE equilibrium populations of the first 
40 l evels o f ortho -H?CO assuming excitation by the 2.7 K background and collisions with 
H 2 . iGreenl (Il99ll ) has calculated excitation rates of these levels for collisions with He taking 
advantage of the sp herical s ymme try of the He potential for kinetic temperatures T = 10 to 



300 K. According to lGreenl (119911 ). excitation rates by H2 collisions could be 2.2 times higher 
than those by He because of the smaller reduced mass and differences in the interaction 
potentials. We will show the results of the calculations under the assumption that the H 2 - 
H2CO collisional rates are t he same as the He-H^CO rates. Probabilities for the radiative 
transitions were taken from Jaruschewski et al. Jl986h . 



The optical depths of the transitions involved in the pumping mechanism are generally 
larger than 1 at high densities and a radiative transfer calculation is required. Two limiting 
approximations in the radiation transport are often considered in molecular clouds: th e 
large velocity gradient (LVG) model and the microturbulent model (ILeung fc Lisztl Il976l ). 
The LVG model assumes that the line profile is dominated by systematic motion of the gas 
while the microturbulent model assumes that the turbulent velocity is much larger than any 
systematic motion. S 140 is likely to have several velocity components, but the existence 
of large systematic motions in molecular clouds and the validity of the LVG model has not 
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been yell establishe d in ot her molecular clouds (e.g. lEvans et al.l (119751 ). IZuckerman fc Evans 
( 11974] ). IZhou et al.l (119901 )). We therefore use the microturbulent model, and for simplicity 
we will use the escape probability formalism to account for photon trapping in a turbulent 
medium in plane-parallel slab of mean total optical depth r t that is perpendicular to the line 
of sight. In this case the A-value of a transition in the equations of statistical equilibrium is 
multiplied by a "loss probability" P(t, r t ) that depends on the mean optical depth r in the 
slab. There are many different ways to define P(r, r t ), which can differ by several orders of 
magni tude for large optical depths. We will use the form suggested by iHummer fc Storey 
(119921 ) for a uniform medium with no continuum absorption: 



P(T,T t ) = m^T) + K 2 (T t -T)) 



where 



dx (j)(x)E2[r(j)(x)} 



(Al) 



(A2) 



and E2 is t he second expon ential integral function. The function ^(t) can be calculated 
from fits by iHummerl (Il98ll ) for a normalized Doppler profile, <j>{x) = eyj£>{—x 2 )/\fK. 



Equation fl All) can be viewed as the single flight escape probability through either side 
of the slab averaged over the line profile. The probability that a photon of an isotropic 
background reaches optical depth r in this slab is also given by equation (1A1I) . and the 
blackbody continuum is thus attenuated by a factor P(t, r t ). The mean optical thickness is 
given by 



dr 



ni 1 



n u /g u 



dz , 



(A3) 



Aun"" \~ ni/gi 

where ai u is the absorption cross section of the transition and Ai/p is the Doppler width. 
We assumed a Doppler width of 2km s" 1 (FWHM = 2Vhi2Av D = 3.3km s" 1 ). The level 
population and its statistical weight are given by n and g respectively with subindex u for 
the upper and / for the lower level. 

The emergent line brightness temperature with subtracted background T is found by 
direct integration of the source function throughout the slab as 



+ 



71 -To 
T [exp( 

hv rT 
~k 



cxp( 



-1] 

-T + r) 



ni/gi 
n u /g u 



-1 



dr 



(A4) 



where 
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m l 



n u /g h 



dz 



, Au D "" V ni/gi 

is the total mean optical depth of a slab of thickness L and T 
total optical depth. 



(A5) 

T t (j)(0) is the line-center 



Equation IA4I gives the correct asymptotic limits: T& — > T when the density goes to 
and Tf, — > T for high densities. The population densities and the optical depths as a 
function of position in the slab in equations (IA1I) and (IA3I) were calculated iteratively. A 
convergence of 10~ 3 K in Tj, was achieved after a few iterations for T < 40 K and H 2 
densities n(H2) < 10 6 cm~ 3 . For higher densities and temperatures the procedure becomes 
unstable. Fig. [5] shows the calculated AT& for a 1 pc-thick slab as a fu nction of the H2 
densit y and a const a nt Ho CO abundance of 2 x 1CT 9 with respect to H 2 (IHasegawa et al. 



(119921 ). iLeung et al.l (ll984f )). The assumed thickness of the slab has an important effect in 
the anomalous absorption as shown in Fig. As the thickness of the slab decreases, the 
effectiveness of the pumping mechanism that cools the line decreases. 
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Fig. 5. — (a) Brightness temperature minus background continuum vs. density for different 
kinetic temperatures T of a 1-pc thick slab, (b) Contour intervals are 0.2 K from AT), = — 1 
to 1 K. (c) Brightness temperature minus background continuum vs. density for different 
slab thicknesses at T = 40 K. 



Garrison et al.l (119751 ) identified some transitions, like In — > 3i2, that produce selection 
effects in the excitation of some levels that cool the H2CO doublets. We have tested our model 
for possible variations of the collision rates. An overall increase of collisional rates by a factor 
of 2.2 decreases AT b in figure by 0.3 to 0.5 K for T > 20 K and n(H 2 ) > 1.6 x 10 4 cm" 3 . 
For lower T and n(H2) there is very little variation in the predicted 
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A. 2. A PDR model for CO 

The UV field has a higher influence on the CO brightness temperature than on the 
H 2 CO brightness temperature. Both molecules are quickly photodissociated near the edge 
of clouds, but the larger abundance of the CO molecule makes its chemistry and interaction 
with radiation more complex. In order to take into account the variation of CO abundance 
along the line of sight, we used the Meudon PDR code to c alculate the CO bright ness 
temperature of a plane-parallel slab irradiated by a UV field (ILe Bourlot et al. Jl993|). A 



detailed description of a revised version of the code is given bv ILe Petit et all (bood ) . 



The sharp separation between the molecular, atomic and i onized emissions suggests 



that the L1204/S140 interface is a PDR viewed nearly edge-on (lHayashi fc Muratalll992l ) 
irradiated by HD 211880. The angle of incidence of the star's radiation on the PDR boundary 
is an unknown parameter but appears to be more-or-less perpendicular. Furthermore, the 
infrar ed embedded sou rces are probably young stars that may also enhance the radiation 



field flEvans et al.lll989h 



We ran the code with its parameters set to represent a plane-parallel slab irr adiated 



from o ne side by a UV field with an enhancement factor x — 200 with respect to the iDraine 



(119781 ) average interstellar radiation field. In a PDR, the gas is heated by the photoelec- 
tric emission from grains and PAH's, H 2 formation in grains, UV pumping in the Lyman 
and Werner bands, gas-grain collisions, photoionization, and photodissociation. As the UV 
radiation is absorbed deeper into the cloud, other processes like cosmic rays and chemical 
reaction energies become important. Cooling is produced by fine-structure and molecular 
line emission. Shocks and turbulence can keep a PDR away from isobaric equilibrium. How- 
ever in our model the temperature and density of the PDR were kept constant in order to 
compare the results with our H2CO model in Fig. [6j We used the chemical network given 
for S140 by the Meudon group at its Internet sitqj, which does not include H 2 CO. Thus the 
CO and H 2 CO calculations represent different models, and Fig. [6] is given only as indication 
of the local conditions that produce the emission and absorption for each molecule. 



A. 3. Photodissociation of H 2 CO 



The n 2 CO molecule is quickly photodissociated into CO and H 2 or H in the average UV 



interstellar radiati on field with a rate of 1.0 x 10 9 exp(— 1.7Av) s 1 (Ivan Dishoeck! Il988l ). 



Keene et al.l (119851 ) estimated that far-ultraviolet radiation (FUV) from the star HD211880 



2 http: / /aristote. obspm.fr/MIS/pdr/exe. html 
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l2 CO H 2 CO 




log n(H 2 cm-') log n(H 2 cm-) 

Fig. 6. — Brightness temperature minus background continuum vs. density for different 
kinetic temperatures T for both the 12 CO(1-0) line and the 6 cm H 2 CO line (Note that 
AT b (k) is plotted "negative" compared to Figure [5]). 



will have an e nhancement factor of Go — 150 wit h respect to the ave rage interstellar field 
( iHabingl Il968l ) at the ionization front, although ISpaans et al.l (119971 ) found that a more 
intense radiation field may be neede d to explain the H2 rota tional emission. H2CO has a 
dissociation energy of 3.61 ± 0.03 eV (ISutoHWang fc Lee 19861). an d a photodissociation rate 
of 1.0 x 10~ 9 sec -1 in the interstellar field ( van Dishoeckj 1988 ) while for CO the values 
are 11.2 eV and 2.0 x 10~ 10 sec -1 respectively. Thus it is possible that the H2CO will be 
selectively photodissociated near the S140 ionization front and the brigh t PDR region, where 
the CO emission peaks. Detailed PDR model calculations by iLi et al.l (120021 ) show that at 
20' from the ionization front, where we observe the H2CO maximum, Go < 20 and Ay ~ 15. 

We added 62 reactions involving H 2 CO and H 2 CO + taken from the UMIST data base 
(Woodall et al., 2007) @ to the chemical network of the Meudon group mentioned above and 
ran a PDR model with Go = 200, T = 40K and constant density of 10 3 cm~ 3 . We found 
that H2CO has significant abundance only at depths of A v > 7 while CO becomes important 
at A v > 4, which shows that photodestruction could explain the offset between the CO and 
the H 2 CO peaks. 



3 http: / / www.udfa.net/ 
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Profile moments in detail 
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Table 2. Profile moments for each position of Fig. [2j The upper limits correspond to 3cr. 
The symbol ". . . " indicates no data were available, while "nf" indicates that data were 
available but no reliable fit could be made. Values marked with "yes" in column 6 are 
plotted as squares in Fig. IU and were used for the least square fits. 



Offset (l,b) 


1UUUX L{ri2^j(JJ 


/T/\ 


1(CUJ 




Included in 


arcmin 


K km s _1 


km s _1 


K km s _1 


km s _1 


nts 


40, 


-70 


< -24.3 


nf 








30, 


-70 


< -18.9 


nf 








20, 


-70 


< -25.8 


nf 








10, 


-70 


< -23.4 


nf 


4.1 ± 0.7 


-7.6 ± 1.3 




o, 


-70 


-42.1 ± 8.4 


-7.5 ± 1.5 


3.5 ± 0.6 


-8.7 ± 1.5 


ves 


-10, 


-70 


-30.3 ± 7.1 


-6.2 ± 1.5 


8.6 ± 0.8 


-8.3 ± 0.8 


ves 


-20, 


-70 


< -21.6 


nf 


< 2.1 


nf 




40, 


-60 


< -22.5 


nf 








30, 


-60 


< -24.0 


nf 








20, 


-60 


< -26.7 


nf 


5.9 ± 0.8 


-7.0 ± 1.1 




10, 


-60 


< -23.7 


nf 


5.1 ± 0.6 


-8.3 ± 1.0 




o, 


-60 


< -19.2 


nf 


< 1.8 


nf 




-10, 


-60 


-27.9± 7.1 


-6.7 ± 0.5 


< 2.1 


nf 




-20, 


-60 


< -28.8 


nf 


< 1.5 


nf 




40, 


-50 


< -25.5 


nf 








30, 


-50 


-109.3 ± 8.9 


-2.6 ± 0.2 


14.8 ± 1.1 


-1.8 ± 0.2 


yes 


20, 


-50 


< -25.2 


nf 


10.5 ± 0.8 


-5.0 ± 0.4 




10, 


-50 


< -21.6 


nf 


8.2 ± 0.6 


-7.0 ± 0.5 




o, 


-50 


< -18.9 


nf 


< 2.4 


nf 




-10, 


-50 


< -16.5 


nf 


< 2.4 


nf 




-20, 


-50 


< -24.6 


nf 


< 4.8 


nf 




40, 


-40 


< -23.1 


nf 








30, 


-40 


< -22.8 


nf 


3.1 ± 0.6 


-6.2 ± 1.2 




20, 


-40 


< -22.2 


nf 


7.2 ± 0.8 


-7.2 ± 0.9 




10, 


-40 


< -23.1 


nf 


4.5 ± 0.5 


-8.8 ± 1.1 




0. 


-40 


< -22.5 


nf 


< 2.7 


nf 




-10, 


-40 


< -26.1 


nf 


< 2.7 


nf 




-20, 


-40 


< -18.6 


nf 


< 3.0 


nf 




40, 


-30 


< -22.5 


nf 


3.4 ± 0.5 


-9.0 ± 1.5 




30, 


-30 


< -22.8 


nf 


4.0 ± 0.7 


-8.9 ± 1.5 
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Table 2 — Continued 



Unset {1,0) 


lOOOx l(H2UUj 




I(CO) 




Included in 


arcmin 


K km s 


km s -1 


K km s 


km s -1 


tits 


20, -30 


< -25.5 


nf 


< 1.8 


nf 




10, -30 


< -30.9 


nf 


5.2 ± 0.6 


-8.6 ± 1.0 




0, -30 


< -20.1 


nf 


< 4.2 


nf 




-10, -30 


< -26.7 


nf 


< 1.8 


nf 




-20, -30 


< -22.2 


nf 


< 2.1 


nf 




40, -20 






< 2.7 


nf 




30, -20 


< -30.9 


nf 


3.5 ± 0.7 


-9.4 ± 1.9 




20, -20 


< -24.0 


nf 


< 2.7 


nf 




10, -20 


-31.6 ± 6.0 


-8.2 ± 1.6 


8.2 ± 1.1 


-8.9 ± 1.2 


ves 


0, -20 


< -29.1 


nf 


8.5 ± 1.4 


-9.5 ± 1.7 




-10, -20 


< -19.5 


nf 


6.8 ± 0.3 


-8.7 ± 0.9 




-20, -20 


< -22.2 


nf 


< 1.5 


nf 




40, -10 






< 1.8 


nf 




30, -10 


< -35.7 


nf 


< 2.4 


nf 




20, -10 


< -16.8 


nf 


< 3.0 


nf 




10, -10 


-43.1 ± 8.2 


-8.1 ± 1.6 


4.7 ± 1.1 


-10.3 ± 2.8 


ves 


0, -10 


-76.2 ± 13.8 


-9.5 ± 1.8 


11.7 ± 0.8 


-9.0 ± 0.6 


ves 


-10, -10 


-82.7 ± 7.7 


-8.1 ± 0.8 


18.6 ± 0.6 


-8.2 ± 0.3 


ves 


-20, -10 


< -24.3 


nf 


< 2.7 


nf 




40, 


< -33.0 


nf 


2.8 ± 0.5 


-11.6 ± 2.3 




30, 


< -28.8 


nf 


< 3.3 


nf 




20, 


-63.8 ± 7.8 


-7.8 ± 0.9 


3.8 ± 1.1 


-8.9 ± 3.0 


yes 


10, 


-93.7 ± 11.3 


-8.7 ± 1.1 


9.4 ± 0.7 


-9.9 ± 0.8 


yes 


0, 


-272.6 ± 9.8 


-7.8 ± 0.3 


21.6 ± 0.7 


-8.0 ± 0.3 


yes 


-10, 


-159.8 ± 8.6 


-8.0 ± 0.4 


46.5 ± 1.6 


-7.8 ± 0.3 


yes 


-20, 


< -24.0 


nf 


10.2 ± 0.7 


-7.7 ± 0.6 




40,+10 






< 1.7 


nf 




30, +10 


< -28.8 


nf 


< 0.9 


nf 




20, +10 


-52.0 ± 8.2 


-8.7 ± 1.4 


5.3 ± 0.4 


-9.1 ± 0.7 


yes 


10, +10 


-114.1 ± 9.7 


-6.2 ± 0.6 


11.1 ± 1.3 


-7.1 ± 0.9 


yes 
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Table 2 — Continued 



Offset (l,b) 


iuuua i^ri2VjWj 


\ v I 




\ v I 


iiici ucicci in 


arcmin 


J.V JA.111 O 


L"ryi c ^ 
ivlll o 


K km s -1 




fits 

11 to 






-70.3 ± 7.7 


-10.3 ± 1.2 


13.5 ± 1.1 


-10.6 ± 1.0 


yes 


o, 


+10 


-247.0 ± fO.f 


-7.6 ± 0.3 


22.8 ± 1.6 


-8.5 ± 0.6 


yes 


-io, 


+10 


-118.9 ± 9.2 


-8.1 ± 0.6 


22.3 ± 2.0 


-8.6 ± 0.8 


yes 


-20, 


+10 


< -19.2 


nf 


< 0.9 


nf 




40, 


+20 






< 1.2 


nf 




30, 


+20 


< -25.2 


nf 


3.1 ± 0.4 


-8.3 ± 1.7 




20, 


+20 


-42.3 ± 9.3 


-7.1 ± 1.6 


6.7 ± 0.7 


-9.9 ± 1.0 


yes 






-82.0 ± 7.4 


-11.3 ± 1.1 


nf 


nf 




10, 


+20 


-49.8 ± 6.9 


-6.9 ± 1.0 


7.0 ± 1.2 


-8.0 ± 1.4 


yes 






-63.1 ± 5.5 


-10.8 ± 1.0 


8.3 ± 0.9 


-11.0 ± 1.4 


yes 


o, 


+20 


-68.3 ± 9.7 


-8.9 ± 1.3 


6.1 ± 0.7 


-10.8 ± 1.3 


yes 


-io, 


+20 


-56.4 ± 6.8 


-8.6 ± 1.1 


5.5 ± 0.5 


-10.8 ± 1.1 


yes 


-20, 


+20 


< -29.7 


nf 


< 2.7 


nf 




40, 


+30 






< 1.6 


nf 




30, 


+30 


< -27.9 


nf 


< 1.6 


nf 




20, 


+30 


< -19.5 


nf 


3.4 ± 0.8 


-11.4 ± 2.7 




10, 


+30 


-54.1 ± 8.8 


-9.3 ± 1.6 


5.6 ± 0.7 


-10.3 ± 1.2 


yes 


o, 


+30 


< -18.0 


nf 


5.0 ± 0.6 


-11.1 ± 1.3 




-io, 


+30 


< -27.0 


nf 


2.0 ± 0.4 


-10.2 ± 2.2 




-20, 


+30 


< -23.4 


nf 


< 3.0 


nf 
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